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In the theory of quantum spin liquids, gauge fluctuations are emergent excitations 
at low energy. The gauge magnetic field is proportional to the scalar spin chirality 
Si ■ S2 >^ S3. It is therefore highly desirable to measure the fluctuation spectrum of 
the scalar spin chirality. We show that in the Kagome lattice with a Dzyaloshinskii- 
Moriya term, the fluctuation in Sz which is readily measured by neutron scattering 
contains a piece which is proportional to the chirality fluctuation. 

It has long been suspected that the spin 1/2 antiferro magnetic Heisenberg model on 
the Kagome lattice may support a spin liquid ground state, i.e., a singlet ground state 
which has no Neel order due to quantum fluctuations. Several years ago the compound 
ZnCu3(OH)6Cl2 (called HerbertsmithiteVwhere the Cu^"*" ions form S = 1/2 local moments 
on a Kagome lattice was synthesized, [s], Q] Despite an exchange constant J estimated to be 
~ 200 K, no magnetic order was detected down to 30 mK. Recent neutron scattering shows 
that the spin excitations are gapless and form a broad continuum. [5[] Thus Herbertsmithite 
has emerged as a strong candidate for the spin liquid state. However, for reasons described 
below, much remains unknown about this material and the connection with theory is tenuous 
at best. There is thus a strong need for more experimental probes to help establish the nature 
of this state of matter. 

Theoretically it has been proposed by Ran et al^ based on projected fermionic wave- 
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functions that the ground state is a U(l) spin hquid, with spinous which exhibit a gapless 
Dirac spectrum. On the other hand, recent DMRG calculations on finite size cylinders show 
strong evidence that the ground state is a Z2 spin liquid, with a substantial triplet gap.Q 
However, the nearest-neighbor Heisenberg model appears to be a very delicate point, be- 
cause a small ferromagnetic next-nearest neighbor exchange J2 ~ — O.OIJ is sufficient to 
destabilize the Z2 state. Isl Meanwhile, more detailed projected wavefunction calculations 
show that the Dirac state is surprisingly stable. Furthermore, the application of a couple of 
Lanczos steps produces an energy quite competitive with the energy of the Z2 state obtained 
by DMRG.[9| Thus while there is general agreement that this ground state is a spin liquid, 
the precise nature of the spin liquid remains somewhat unsettled. 

Experimentally it is known that about 15% of the Zn {S = 0) ions which are located 
between the Kagome planes are replaced by 5* = 1/2 Cu ions. It has been argued that there 
is not much Zn substitution for Cu in the Kagome planes, jl(| so that the disturbance of the 
Kagome structure may be minimal. However, much of the low energy excitations measured 
by thermodynamic probes such as specific heat and spin susceptibility are dominated by the 
local moments between planes. Furthermore, due to spin orbit coupling, we expect deviation 
from the Heisenberg model. To first order in the spin orbit coupling constant A, we expect 
Dzyaloshinskii-Moriya (DM) terms of the form 

"^DM = Dij ■ Si X Sj (1) 

where the DM vector Dij is located on bond (r/). Since Dij = —Dji, the Dij vectors depend 
on the convention of the bond orientation. 11[ For a given convention the DM vectors are 
shown in Fig. 1. There is also an exchange anisotropy term Jzz 7^ Jxx- However, the 
difference is second order in A and may be less important than the DM term. The out-of- 
plane DM term (-D^) has been estimated to be about 8% of J. Due to the delicate nature 
of the ground state of the Heisenberg model explained above, it is not at all clear that the 
nearest-neighbor Heisenberg result applies to the Herbertsmithite. 

The defining character of a quantum spin liquid is the emergence of exotic particles such 
as spinous which carry S" = 1/2 and the associated gauge fields. [12|] In the U(l) spin liquid, 
the gauge field is gapless whereas in the Z2 spin liquid the gauge field is gapped. The gauge 
field is defined by the phase aij of the spinon hopping matrix element te*°'^ on link ij. It is 
a compact gauge field and the spin liquid corresponds to the deconfined phase of the gauge 
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FIG. 1: The DM vectors on the Kagome lattice for Herbertsmithite. The arrow specifies the order 
of the operator Si x Sj. (adapted from Ref. 11) D is the z component while D' is the in-plane 
component. 

field, so that the compactness may be ignored in the long wavelength limit and aij may 
be replaced by a continuum field a{r). The gauge invariant quantities are the gauge field 
b = V X a which is in the z direction and the gauge electric field e = — Voq + ^ which lies 



in the plane. The physical meaning of the magnetic flux has been extensively discussed. |13[] 
If $ is the flux through a plaquette, sin$ is one half of the solid angle subtended by the 
spins along the plaquette. For a three site triangle, we have 



sin$ = ^Si -(82x83) 



(2) 



The quantity Si- S2X is known as the scalar spin chirality. Thus the fluctuation spectrum 
of the gauge magnetic field is proportional to the fluctuation of the spin chirality. It is highly 
desirable to measure this correlation function because it gives information on the gauge 
fluctuation and can help to distinguish different spin liquids. However the measurement of 
a correlation of a product of 3 spin operators is a daunting task. A method to measure 
the vector chirality Si x ^2 has been suggested by Maleev (ref 14), but that does not 
apply to scalar chirality. Shastry and Shraiman (ref. 15) have suggested measuring chirality 
fluctuations using Raman scattering, but that contains information only for very small q. A 
proposal to measure this using resonant X-ray scattering has been made. IGj However, the 
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energy resolution of this technique is currently limited to 20 meV or so, which is on the scale 
of J. 

It turns out that in the Kagome lattice we can turn the DM term to our advantage 
and achieve a simpler measurement of the chirality correlation. The hint comes from a 



recent paper by Savary and Balents, 171] who pointed out that in a certain pyrochore spin- 



ice material, the gauge field (in their case the electric field[l8|]) is proportional to and 
its fluctuation can be directly measured by neutron scattering. This system is treated in 
the strong spin-orbit coupling where J rather than spin is a good quantum number, but 
this work raises the possibility of finding the same proportionality in the presence of weak 
spin-orbit coupling. We find that it is indeed the case for the Kagome lattice with the DM 
vectors shown in Fig. 1. In the Kagome lattice each site is connected to two triangles. Let 
us consider the site labeled 1. The total chirality through the two attached triangles is 

Xi = 51 • 52 X 53 + 5i • 54 X 55 (3) 

However, due to the DM term, 

(52 X 53) = aD2'i 

(54 X 55) = aD^^ = -aD^A = aD^^ (4) 

where a is a constant and the last step is by inspection from Fig. 1. For I-D23I ^ J ■, "we 
expect a ~ J. The important point is the contributions from the two triangles add and we 
find a linear coupling between xi ^"^^ 5i ■ -D23- If we average over corners of the triangle 123 
and define x = |(Xi+X2+X3)) it is clear that the in-plane component of the D vectors cancel 
and a local fluctuation in Sz{r, t) induces a local fluctuation in x(r, t) with a proportionality 
constant of {2aDl^). This suggests that a measurement of the {Sz{r,t)Sz(S^,Q)) correlation 
function will contain a piece which is proportional to the chirality correlation (x(^) ^)x(O))- 
A more formal argument proceeds as follows. Let \a) and E{a) denote the exact eigenstates 
and energies of the Hamiltonian Hq without the DM term and we perturb in H-qm- We 
are interested in a subset of the excited states, denoted by |q;^) which are connected to the 
ground state |0) by the operator xii i-6-, («xlxi|0) 7^ 0. On the other hand, the operator 
5*2 (r) has no matrix element to these states, i.e., {a^\Sz{r)\Q) = 0. The gauge magnetic 
field spectral function is proportional to the chirality spectral function 

S^iq,u) = J2\Mx,m'6iuj - {E{a^) - E(0)) (5) 



5 



Now we turn on the Dirac terms. The neutron scattering cross-section is proportional to 

S{q,u) = J2 \{f\S.{qm'S{uj - {Ef - E,)) (6) 
/ 

where |/), E{f) are the exact eigenstate and energy of the total Hamiltonian. We now 
compute the correction to the matrix element {f\Sz\i) to first order in H-dm for the state 
1/;^) which derives from the state i.e., 

l/x) = K)+}:^^( (7) 



(8) 



and similarly for Since the zeroth order term vanishes, we find 

/f \Q ( \r\ \ {(^x\^GM\<y){<y\Sz\0) , {a^\S^\a){a\HDM\0) 

^^'^^''^"t^L E{a,)-Eia) + E(0) - 

Since |0) and \a^) are total spin singlets, we argue that by choosing a different spin 
quantization axis, Sz can rotated to Sx in Eq.(8) and |a) are spin triplet states. If the 
spin triplets have a large gap we may replace £'(0) — £'(0;) in the second term by 
— Af. If we are interested in low energy modes of the chirality fiuctuations such that 
u = E{a^) — E{0) ^ E{a) — E{0) ^ At, we may likewise replace that energy denomi- 
nator in the first term by — Aj. Then the sum over |a) can be done and it is clear that 

9 n 

= -J2^MS.ir,)z-S{rj) X 5(n)|0) (9) 

(Jk) * 

where the sum is over sites j, k which are connected to site 1 as corners of a triangle, 
i.e., {j,k) = (2,3) or (4,5). Putting this matrix element into Eq.(6) we see that neutron 
scattering indeed has a piece which couples to chirality fiuctuations, with a form factor given 
t'Y I S{jfc) least in the limit where the triplet gap is large. We expect this coupling 

will continue to hold even if the assumption of a large At fails, as in the Dirac spin liquid 
case, but the form factor may acquire some quantitative differences and perhaps some q 
and u dependence. This will be the case if the contribution to the sum in Eq(8) from low 
energy (near gapless) triplet excitations are reduced due to restricted phase space and we 
can approximately replace the energy denominator by an average triplet excitation energy. 

As au example, .f the ground state is deser.bed by the U(l) D,ra. spin hqu.d^he eorr. 
lation of the gauge magnetic field is expected to be (in the RPA approximation pj,]) 

(10) 
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where 9 is the step function and v is spinon velocity. The neutron scattering intensity which 
couples to the Sz channel is expected to have a piece given by Eq.(lO), with the intensity 
reduced by (2aDfj)^. If the ground state is a Z2 gapped spin liquid, the gauge fluctuation 
is gapped and one may expect to see a gapped mode instead of Eq.(lO). As pointed out 
by Savary and Balents,jl7] the dependence in Eq.(lO) implies that the gauge fluctuation 
structure factor vanishes as g — )■ 0, in contrast with the spin wave structure factor which 
diverges as 1/q in this limit. While this helps identify the gauge fluctuation, it also makes 
its detection more difficult. 




FIG. 2: The distortion of the Ir (solid) and oxygen (open) ions in the Ir-0 plane of Sr2lr04. 
(adapted from Ref. 20) 



The linear relationship between chirality fluctuation and Sz fluctuation is rather special 
to the Kagome lattice. Such a coupling does not exist for the square lattice, for instance, 
for slowly varying chirality fluctuations. Consider magnetic ions on a square lattice each 
surrounded by corner sharing oxygen cages which rotate about the z axis in a staggered 
manner. (Sr2lr04 is such an example [20! shown in Fig. 2.) In this case each magnetic ion is 
connected to 8 triangles whose opposite side is a nearest-neighbor bond with a DM vector in 
the z direction. The gauge magnetic flux through each square is given in terms of the scalar 
chirality of the triangles which split up the square, [isl It is easy to see that the contributions 
from the two triangles 123 and 134 which split up a square cancel each other and there is 
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no linear coupling between the gauge flux and in this case. 

In conclusion, we have shown that in the Kagome lattice, the DM term leads to a linear 
coupling between the Sz fluctuation with the spin chirality fluctuation. It will be interesting 
to see if neutron scattering can give information on the chirality fluctuation and shed light 
on the nature of the spin liquid in Herbertsmithite. 
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